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The Hahm-Kulsrud (HK) [T. S. Hahm and R. M. Kulsrud, Phys. Fluids 28, 2412 (1985)] solu- 
tions for a magnetically sheared plasma slab driven by a resonant periodic boundary perturbation 
illustrate fully shielded (current sheet) and fully reconnected (magnetic island) responses. On the 
global scale, reconnection involves solving a magnetohydrodynamic (MHD) equilibrium problem. In 
systems with a continuous symmetry such MHD equilibria are typically found by solving the Grad- 
Shafranov equation, and in slab geometry the elliptic operator in this equation is the 2-D Laplacian. 
Thus, assuming appropriate pressure and poloidal current profiles, a conformal mapping method 
can be used to transform one solution into another with different boundary conditions, giving a 
continuous sequence of solutions in the form of partially reconnected magnetic islands (plasmoids) 
separated by Syrovatsky current sheets. The two HK solutions appear as special cases. 
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I. INTRODUCTION 

Recently, there has been renewed interest in the sec- 
ondary tearing instability of high-Lundquist-number cur- 
rent sheets [JJ, called the "plasmoid instability" [2J. Nu- 
merical simulations, supported by heuristic scaling ar- 
guments [5] , demonstrate that if the Lundquist number 
(S) based on the length of the current sheet exceeds a 
threshold [H [5] , the instability breaks up a Sweet-Parker 
current layer into a sequence of magnetic islands sepa- 
rated by segments of current sheets, and evolves into a 
new nonlinear regime of reconnection in which the recon- 
nection rate becomes nearly independent of S EHS]. 
These simulation results suggest that there might exist 
partially reconnected plasmoid solutions of the magne- 
tostatic equilibrium equations in which plasmoids exist, 
separated by segments of current sheets. In this pa- 
per, we show that such solutions can indeed be con- 
structed within the framework of the Hahm-Kulsrud- 
Taylor (HKT) model, described below. 

The HKT model, developed by Hahm and Kulsrud 
[TU] following a suggestion by J. B. Taylor, considers the 
response of a plasma slab with a sheared unperturbed 
magnetic field to a resonant perturbation applied at the 
boundaries x — ±a. In Cartesian coordinates x, y, Z 
the magnetic field B is represented as VzxV^i,)/) + 
B z (x,y)Vz. The unperturbed "poloidal" flux function 
vb is ipo(x) = B a x 2 /2a, where the constant B a is the 



strength of the poloidal magnetic field at the boundary 
x = a. 

In Ref. [10] the "toroidal" field B z was assumed to be 
effectively constant and much larger than B a to allow 
incompressibility to be assumed during a discussion of 
reconnection dynamics, but in this paper we will be con- 
cerned only with finding static equilibrium solutions so 
this assumption is not necessary. Instead, we regard the 
B z profile function F(vj) in the Grad-Shafranov equa- 
tion as free to choose, and assume it is chosen so that 
B 2 /2^jLq+p is linear in ijj, where p(ij>) is the pressure (and, 
for SI units, (io is the permeability of free space). This 
makes the Grad-Shafranov equation linear (though inho- 
mogeneous) allowing analytic solutions to be obtained. 




* robert.dewar@anu.edu.au 
t abhattac@pppl.gov 
T rkulsrud@pppl.gov 



FIG. 1. Contours of tp — tpo + aip for the standard illustrative 
values used throughout this paper: a = 0.1, units such that 
a — B a = 1, and wavelength A chosen to equal a, so k = 2-7T. 
This illustrates the effective rippling of the upper boundary 
(axes scales equal). Cuts corresponding to the case L = A/6 
are also shown as thick gray lines on the j/-axis. (Color online.) 

The perturbed flux function at the boundaries is as- 



sumed to be 




ip(±a,y) 


— \aB a + B a 8 cos ky 




= aB a (| + a cos ky) 



(1) 



where a = S/a is a dimensionless parameter measuring 
the strength of the perturbation, 8 being the amplitude of 
a notional small boundary ripple of wavelength A = 27r/fc, 
the given planar boundary conditions being, to linear or- 
der in a < 1, equivalent to a symmetric geometric rip- 
pling of perfectly conducting bounding walls. 

This rippling is illustrated in Fig.nl), though the value 
of a = 0.1 we use in our standard illustrative case is 
rather too large for the ripples to be even approximately 
sinusoidal. However, in this paper we are do not really 
need to ripple the boundary and take the HKT bound- 
ary conditions Eq. (IT]) as exact, so that a need not be 
infinitesimal. Thus, because of the assumed linearity of 
the Grad-Shafranov equation mentioned above, we have 
exact linearity and can write 



full current sheet is one of a continuous infinity of full 
current sheet solutions differing by the strength of a con- 
stant inten sity of current in the sheet, the HK solution 



ip = ipQ(x) + atp(x,y) 



(2) 



where tp is independent of a and obeys the boundary 
conditions ip(±a,y) = aB a cos ky. 

Ideal magnetohydrodynamics (MHD) cannot predict 
the timescale for magnetic reconnection (magnetic field 
changes that violate the topological frozen-in-flux con- 
dition). However, on a long enough length scale and a 
short enough time scale, the intermediate states in the 
evolution of driven [TH - TH] or spontaneous [T5] magnetic 
reconnection can be described as a continuous sequence 
of "global" MHD equilibrium states (i.e. states that sat- 
isfy the boundary conditions and internal force balance) . 
In the present paper we do not seek to describe the re- 
connection process in detail, merely to find analytically 
solvable Grad-Shafranov equilibria that plausibly illus- 
trate a possible reconnection scenario. 

In Sec. [IT] it is pointed out that the Grad-Shafranov 
equation in slab geometry can include current sheets 
in two ways — either as a superposition of (^-function 
current-density sources or as cuts in the x, y plane, the 
latter being the viewpoint used in this paper. Force 
balance provides the boundary conditions on the cuts. 
Current profiles are given such that the Grad-Shafranov 
equation in slab geometry becomes a linear Poisson equa- 
tion and our definition of the HKT equilibrium problem 
is made precise. 

Hahm and Kulsrud [10] found two exact MHD equilib- 
rium solutions, one involving a full current sheet cover- 
ing the plane x = and one describing a magnetic island 
with no current sheet. The current sheet solution may 
be viewed as representing how a shielding current |16[ 
c.g] initially arises in order to prevent reconnection af- 
ter a resonant perturbation is turned on, the magnetic 
island solution being interpreted as the end state after a 
sufficient time has elapsed that reconnection has run its 
course and an island has "opened." In Sec. |IH| we review 
the Hahm-Kulsrud (HK) solutions and show that their 



(■0i in Sec. Ill) being the one with zero net current. 
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FIG. 2. Contours of — ipeut, for the partially-reconnected 
plasmoid case L = A/6 with 7s = 1. Here cu t = —0.000436 
is the value of ip on the cuts shown as black lines on the y-axis. 
(Color online.) 

The existence of solutions to the HKT equilibrium 
problem with either a full current sheet or no current 
sheet raises the possibility that there may be more so- 
lutions, intermediate between the two solutions found in 
Ref.[10l As discussed above, current sheets can be unsta- 
ble to the formation of plasmoids embedded in the cur- 
rent sheet. Thus we seek ideal equilibria that illustrate 
topologically a scenario for the decay of the shielding 
HKT current sheet via a plasmoid mechanism of island 
formation, continuously connecting the two extreme HK 
solutions. While our solution is not exact, except for 
the two HK limiting cases, it satisfies Eq. (IT]) very ac- 
curately [error 0(e~ 3ka )] so could presumably be made 
exact by a small perturbation of our ansatz. A typical 
plasmoid case is depicted in Fig. 151 where current sheets, 
the black horizontal lines of length 2L = A/3, alternate 
with plasmoids of width 2A/3. A magnified view of a 
typical current sheet end-point for the case L = A/2 is 
shown in Fig. [3] 
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FIG. 3. A magnified view of the ip contours in the vicinity of 
the junction of a current sheet and a plasmoid, for the case 
L = A/4, 7s = 1, showing the Y-point structure obtained for 
this value of 7s (axes scales equal). 

In Sec. |IV| we review Syrovatsky's complex variable ap- 
proach to finding an analytic solution of Laplace's equa- 
tion that represents a large-scale view of a Sweet-Parker 
current sheet [TJJ e.g.] (see Fig. [9]). Syrovatsky's solu- 
tion was obtained by using a conformal mapping from 




the simpler solution for the field around a neutral point 
in the poloidal field. 

In Sec.fVlwe introduce a new, periodic conformal map- 
ping to transform the Syrovatsky solution into a plasmoid 
solution of the HKT equilibrium problem. 

In Sec. |VI| we analyze the mismatch between the 
boundary condition Eq. (TTj) and ip obtained from our 
conformal mapping ansatz. We present this error both 
graphically, vs. y and L in a typical case, and also give 
an analytic expression for the first nonvanishing term 
(6th order!) in an expansion in the small parameter 
e = exp(— ak/2). Possible further improvements and 
applications of the plasmoid scenario are discussed in 

Sec. Ivnl 



II. GRAD-SHAFRANOV EQUATION WITH 
CURRENT SHEETS 



For analyzing equilibria with ideal (zero thickness) cur- 
rent sheets, the force-balance condition is best written in 
the conservation form 



B 2 ^ 

2p 



I 



BB 

Mo 



= 



(3) 



where p is the plasma pressure and B is the magnetic field 
(SI units). In regions where p and B are differentiable, 
this implies the force balance condition, jxB = Vp, 
where j = ^'(VxB) is the plasma current. However, 
in the neighborhood of an ideal current sheet, p and B 
are not everywhere differentiable and we need to use gen- 
eralized functions, like the Dirac delta function S(s), to 
find weak solutions of Eq. ([3| . 

However, we can avoid using generalized functions ex- 
plicitly by cutting the x, y plane along its intersections 
with current sheets and solving on the cut plane with ap- 
propriate boundary conditions on the cuts. The bound- 
ary conditions on the two sides ± of a current sheet are 
found to be |17l Appendix A] the tangentiality conditions 



n-B- 







and the pressure-balance jump condition, 

B 2 



P- 



2^0 







(4) 



(5) 



The first condition implies that an equilibrium current 
sheet must be a tangential discontinuity in B. 

In the case of a cylindrical or slab plasma of arbitrary 
cross section (independent of z) , a general representation 
for the equilibrium magnetic field is 



B = VzxVip + F(ip)Vz , 



(G) 



where 0(x, Tj) is the flux function defined in Sec. Kj x,y,z 
being Cartesian coordinates with the z-axis in the sym- 
metry direction. The first of Eqs. (pj implies that the two 
sides of an equilibrium current sheet are level surfaces 



of ip± (in fact -0 must be continuous across the current 
sheet, 0_ = ip+ , to avoid infinite poloidal magnetic field 
there) while the second gives 



P- 



|V0| 2 +F 2 



2p 



= 



(7) 



Taking the curl of Eq. wty we find, everywhere except 
on a cut, 



Moj = VVVz + F'()/))V^xV2 



> 2 +d 2 , 



(8) 



and 



where V 2 is the 2-dimensional Laplacian, d 2 
F' = dF/dip. 

Summarizing, the equilibrium condition Eq. ([3]) is sat- 
isfied if and only if the two-dimensional Grad-Shafranov 
equation for axisymmetric static MHD equilibria, 



v 2 v + fy [moXV) + \ F (i>f] = o 



(9) 



is satisfied everywhere except on cuts, where the current 
sheet force-balance condition Eq. ^f\ applies instead. 

Equation ([9]) is in general nonlinear, but consider the 
special case fiop(ip) + |-F(0) 2 = const — (B a /a)ip for 
which Eq. M becomes a Poisson equation, V 2 = B a /a, 
linear in 0. This can be solved as a linear superposition, 
■0o + 0h where 0h is a harmonic function, i.e. a solution 
of the Laplace equation 



V 2 h 







(10) 



determined by the boundary conditions. Comparing with 
Eq. (12| we identify V'h with atp. The wall boundary con- 
ditions Eq. (IT]), current sheets on the y-axis, the assumed 
form of ipoizji an d Eq. (10), make up what we call the 



HKT equilibrium problem, whose scope we expand by 
considering a wider class of current sheet cuts in the do- 
main on which Eq. ( 10 ) is to be solved. This has the 



consequence that tp cannot be assumed necessarily to be 
sinusoidal in y. 

Noting that [0] = we see from the assumption 
£top(0) + \F{' I P) 2 — const — (B a /a)ip that the equilib- 
rium jump condition Eq. ufh simplifies to 



|V0| 2 1 = 



(11) 



Although is continuous, and by Eq. (11), |V0| is con- 
tinuous, V0 can be discontinuous, its jump giving the 
intensity j t of the (5-function component of the current 
j = VxB/^ , 



J = 



1 

/in 



F'(i/>)V0xVz 
B a 



+ ld x il>}5(x))Vz 



(12) 



at each point on the current sheets on the y-axis, whence 
3*{y) = Idxipj/fM)- 



WAJ 



It is easily verified that 
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FIG. 4. Jumps in the gradient of ip across the y-axis, giv- 
ing the intensity of the sheet currents by Eq. (12 1, for the 
original HK fully screened solution 7s = (dashed) and the 
modified HK solution 7s = 1 (solid) for which the current 
sheet intensity never goes negative. 



III. GENERALIZED HK SHIELDING 
SOLUTIONS 



Hahm and Kulsrud [10 found two solutions of 
the HKT equilibrium problem, a shielding current 
sheet solution ipi(x, y) — aB a |sinh(fcx)| cos(fcy)/sinh(fca), 
and a fully developed island solution ipu(x,y) — 
aB a cosh(/cx) cos(fcy)/cosh(£;a) with no current sheet [T5] 
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FIG. 5. Contours of ip — ipcut near the y-axis for the original 
HK fully screened solution ip — ipo + anpx (ip\* for 7s = 0), 
showing islands with separatrices joining points where the 
sheet current changes sign. 
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FIG. 6. Contours of ip — ipcut near the y-axis for the modified 
fully screened solution ip — ipo -+- otipi* with 7s = 1, showing 
the simple topology of the magnetic surfaces in this case. 



ipi,{x,y) = aB a 



|sinh(fea;)| cos(ky) + jsk(\x\ — a) 
sinh(fca) 



(13) 



where 7s is an arbitrary constant, also satisfies the HKT 
equilibrium problem. The inclusion of the 7s term rep- 
resents a small but important generalization of the HK 
shielding current sheet solution that allows the dc level 
of the current in the sheet to be adjusted, as illustrated 
in Fig. [4] 

This has a rather profound influence on the topology 
of the ip contours, as illustrated in Figs. [5] and [6] where 
it is seen that the HK solution ip has saddle points at 
the current-reversal points x = 0, y = A/4 ± nA/2, n = 
0,1,2,..., where the ip = contour bifurcates off the 
y-axis, forming magnetic islands. 

In the following sections we obtain these solutions as 
limiting cases of a new family of solutions ip obtained 
using a conformal mapping approach [19 , which relies on 
the facts that the real or imaginary part of any analytic 
function, .Fg(£), £ = x + iy, is harmonic, and that the 
composition of two analytic functions is itself analytic. 



IV. THE SYROVATSKY CURRENT SHEET 



The double- valued analytic functions [20] 
F s (u) = sgn(Reu) lu{u 2 + 1) 1/2 

+ 7s ln[ U +( M 2 + l) 1 / 2 ]} , 



Fg (u) — sgn (Re u 



2u 2 



•7S 



(u 2 + 1)V2 



(14) 
(15) 



defined on the complex w-plane with two Riemann sheets 
joined by a cut joining branch points u = ±i, may be used 
to define the harmonic functions ips(x',y') = KeFs(x' + 
iy') and d x iips{ x ' >v') — Re-Fg(x' + iy')- The former can 
be interpreted as the flux function for a Sweet-Parker 
current sheet positioned on the cut between y' — ±1 and 
the latter gives the intensity of the current sheet on the 
cut. 

The form in Eq. ( 14 ) is that used in Ref . [5] with square 
root z 1 ' 2 and natural logarithm In z defined as usual on 
the complex z-plane cut along the negative real axis. The 
step function factor, sgn (Re u) = ±1, is needed to make 
the cut a straight line joining y' = ±1 rather than make 
two cuts radiating outward to infinity. The function de- 
fined by Eq. (14 1 is obtained by the following substitu- 
tions in Syrovatsky's [20] Eq. (44), 



n i->- 2, b i->- i, A >-> — ,a n- 2, r n- 2^73,2 n- (' 



It is seen from Eq. ([15]) that in general the sheet cur- 
rent intensity diverges at the endpoints y' = ±1, but the 
choice 7s = 1 [301 Eq. (47)] makes the current go continu- 
ously to zero at the endpoints, giving rise to the Y-points 



seen in Fig. [91 This figure shows contours of tps(x' , y') in 
the y'jx' plane, the thick horizontal line indicating the 
cut/current sheet and the thinner continuous lines the 
magnetic field lines. The ellipses and hyperbolae form a 
visualization of a periodic conformal map shown in the 
next section to convert Syrovatsky's single current sheet 
solution into a plasmoid solution of the HKT equilibrium 
problem. 
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FIG. 7. Jumps in the gradient of ip across the j/-axis for the 
half-screened plasmoid case L = A/4. The dashed curve is for 
the case 7s = 0, showing current reversal within the current 
sheets and singularities at the ends, and the solid curve is for 
the case 7s = 1 in which the current sheet intensity never 
goes negative and joins continuously to the zero values within 
the plasmoid regions. 



SHINUSOIDAL TRANSFORMATION 



The analytic function 



/(C) 



sinh(|fc£) 
sin(|fci) 



(16) 



will be used to map the single current sheet in the 
Syrovatsky solution Eq. (14 1 to a periodic sequence of 
current sheets of length 2L, replacing the X-points at 
y = n\, n = 0, ±1, ±2, . . ., in the fully reconnected (mag- 
netic island) solution of Ref. [TU1 where A is the wave- 
length of the boundary perturbation. 

Its appropriateness for this purpose will be verified be- 
low a posteriori, but as partial motivation for this ansatz 
we note some useful properties of /: 

• f(x) = s'mh(^kx)/ sin(^kL) 
f(iy) = isin(ifcy)/ 'sin(lfcL) 

f(x ± \i\) = ±i cosh(^kx) / sin(^kL) 

• f( x + iy) is periodic in y with wavelength 2A, but 
f(x + iy) 2 has wavelength A. 

• The only zeros of /(£) are at £ = inX, i.e. at the 
X-points of the fully reconnected solution, while at 
the first O-point, y = A/2, f(iX/2) = i/ sm.{\kL) 
ranges from i in the fully reconnected case, L = 
A/2, to ioo in the fully shielded case, L = 0. 




FIG. 8. Contours of ip — ^o+ctip for a plasmoid case L = A/4, 
7s = 1, with magnified scale on a;-axis. Mapped images of the 
vertical (orange) and horizontal (green) mesh lines appear in 
Fig. [9] (color online). 



• The double- valued function sin(|fcL)[l + 
/(y) 2 ] 1/2 = [(sinffcL - sin |fcC)(sin \kL + 
smifcC)] 1 / 2 has, within the strip —A/2 < Im£ < 
A/2, two branch points, located at the endpoints 
( = ±L, of the sought-for current sheet,. 




FIG. 9. Contours (black) in the y' , x' -plane of the Syrovatsky 
solution il>s(x',y') — HeFs(x' ,y') given in Sec. IV in the case 
7s = 1- The orange hyperbolae and green ellipses are the 
images of the rectangular mesh in Fig. [8] extended to the 
range x 6 [—0.2,0.2], under the periodic conformal map x' + 
iy' — f( x + W) defined in Sec. [V] (axes scales equal; color 
online). 



We now use the conformal map £' = /(C) to trans- 
form the Syrovatsky function to a function of £, -Fs(C') — 
i ? S°/(C)) that provides the harmonic function ijj(x,y) 
through the equation 



i>(%, y) = aB a 



ReF s (C) 



(17) 



where C = f( x + w) an d the constants c and d are 
determined by requiring that the boundary condition 
■0(±a, y) = clBq cos ky be satisfied to a good approxi- 
mation (see Sec VI). 

Figure [2] illustrates how this transformation results in 
a typical plasmoid structure for < L < A/2. A graph- 
ical visualization of the transformation may be had by 



comparing the mesh in Fig. [8] with its image in Fig. [9] 
Figure [l] illustrates the decay of the ripple away from the 
boundary before it is amplified by the resonance effect 
near the poloidal field reversal region, \x\ ~ 0, as seen in 
Figs. [5] and [6| 

Figure [7] shows the current sheet intensity for the cases 
case 7s = 1 and case 7s = 0, showing the singular be- 
havior inherited from the Syrovatsky solution Eq. ( 15 ) 
in the latter case. Figure [3] verifies that the case 7s = 1 
leads to the typical Y-point magnetic surface behavior 
expected of a Sweet-Parker current sheet, whereas Sy- 
rovatsky [20l Fig. 3] showed that the case 7s = leads 
to reentrant ^-contours producing a cusp pointing away 
from the current sheet. 

It ma y be shown analytically that the solution 
Eq. (17) reduces to the HK island solution tp(x,y) — 
aB a cosh(/cx) cos(fcy)/cosh(A;a) as L — > 0, while it reduces 
to the generalized HK current sheet solution Eq. |l3| ) as 
L ->■ A/2. 

iA(«,v) - cos&y 



rs=o 
ys=i 



FIG. 10. The residual boundary error ip(a,y) — cosky vs. y 
for the case L — A/4. 



in order to determine the two constants c and d in 
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FIG. 11. The residual boundary error maximized over y, i.e. 
|$(a,0)-l|, vs. L. 



Eq. (17). The two conditions are 



e„(0) - £ Q (A/2) = and e„(0) + e„(A/4) = . (19) 

We can now verify a posteriori that the boundary con- 
ditions are satisfied to high accuracy for all y in typical 
cases. For instance, in Fig. [l0]we plot e a (y) for the same 
case as shown in Fig. [8] and see that the conditions in 
Eq. ( 19 ) null out any constant error and the fundamental, 



cos ky, leaving only a second harmonic error proportional 
to cos 2ky, with an amplitude that is extremely small in 
the case studied. 

Noting that the error is an even function, periodic in 
y, we see that y = is a maximum point of the absolute 
value of the error. In Fig. 11 we plot this maximum 
error vs. the halfwidth, L, of the current sheet. This 
figure shows that the error is zero for the two HKT cases 
L = (complete reconnection) and L = A/2 (complete 
shielding), and nowhere gets much larger than it does in 
the typical intermediate case depicted in Figs. [8| and [T0| 

In the plots Eqs.[l9]have been solved numerically, but 
to understand why the error is so extraordinarily small 
it is instructive to perform a perturbation expansion in 
e = exp(— ak/2) (— 0.0432 ... in the reference case used 
in this paper), which allows us to separate those terms 
that are large at x = a from those that are small. For 
instance, 



f(a + iy) = 



exp(|ifcy) — eexp(— \iky) 
2sin(|fci) 



(20) 



Expanding Eq. (17) in e and solving for c and d such 



that the constant and cos ky terms in e a (y) vanish we 
find, to first nonvanishing order, the residual boundary 
error 

aB n 



3ak sin 2 kL 
< [2 - 7s - (2 - 37s) cos kL] cos 2ky (21) 



ip{a,y) - aB a cosky = —e 



which is 0(e 6 ), thus explaining the smallness of the error 
found numerically. 



VII. CONCLUSION 



VI. BOUNDARY ERROR ANALYSIS 



The current sheet force-balance requirement Eq. (11) 



is ensured by the restriction of the cuts to the y-axis, and 
the assumed symmetry about this axis, but the bound- 
ary conditions Eq. (fTl) are not imposed a priori for all 
y. Instead we impose only two conditions involving the 
boundary error function, 



£a{y) 



ReFs[/(a + iy)] — d cosky 



(18) 



This paper demonstrates that solutions to the global 
ideal-MHD equibrium problem are far from unique when 
interior current sheets are allowed. We have made no at- 
tempt to link the members of our equilibrium sequence by 
applying constraints consistent with almost-ideal-MHD 
time evolution, so this work, of itself, cannot be regarded 
as a study of reconnection. However, it is highly sugges- 
tive that evolution through a plasmoid phase represents a 
topologically reasonable mechanism for an initial shield- 
ing current sheet to open into a magnetic island. To 
establish this scenario as a reconnection mechanism, two 
approaches appear promising, both applying a subset of 
the ideal-MHD constraints: 



1. A maximally constrained or almost-ideal MHD 
approach assuming ideal MHD applies locally 
throughout the evolution, except as a plasma passes 
through a Sweet-Parker current sheet where the 
frozen-in-flux constraint is relaxed and reconnec- 
tion can occur. While respecting the detailed 
physics of the process, it is not amenable to the 
conformal mapping approach we have used to find 
analytical solutions as it does not preserve the con- 
dition of linearity of nop(ip) + F(ip) 2 /2 assumed 
at the beginning of Sec. |VJ Furthermore, it implies 
current sheets on the plasmoid separatrices, so that 
the simple cut structure of the Syrovatsky solution 
does not apply [121 US]. Thus a completely different 
method of analysis would need to be applied. An 
interesting approach has been discussed by Kulsrud 

EH- 

2. A minimally constrained or relaxed MHD approach 
based on a generalization of Taylor [22 relaxation 
to include more ideal-MHD invariants than the 
magnetic helicity constraint assumed by Taylor, 
but only a sufficient number to capture the qual- 
itative essence of the evolution (cf. [531 HI])- A 
noncanonical Hamiltonian approach has recently 
been developed [33] in which the ideal-MHD con- 
straints appear as Casimir invariants. In this work 
it was shown that bifurcation of a cylindrical Tay- 
lor relaxed state to a helical relaxed state can be 
frustrated by introducing a singular Casimir invari- 
ant corresponding to the shielding HKT current 
sheet, the magnetic field everywhere else in the 
plasma being given by the linear Beltrami equa- 
tion, VxB = /xB, found by Taylor. This suggests 
seeking, in slab geometry, a sequence of plasmoid 



solutions analogous to those found in the present 
paper, especially in the limit (i — >• where the Bel- 
trami field reduces to a harmonic field correspond- 
ing to ip. 

Using either approach to generate an equilibrium se- 
quence with fixed boundary conditions, its applicability 
as a physically plausible reconnection scenario could be 
determined, without the necessity of resolving the cur- 
rent sheets into finite-width tearing layers, simply by 
showing that the plasma potential energy W = /[p/(7 — 
1) + B 2 /2[io]dV [2S] decreases monotonically along the 
sequence, the final state being a minimum of W. Pre- 
sumably, if two sequences are parametrized by their re- 
connected fluxes and the graph of the potential energy of 
one lies below that of the other, then the first sequence 
is physically preferred. This could be used to determine 
when and if the symmetry-breaking plasmoid evolution 
found in Ref. [5] can occur, rather than the symmetric 
evolution assumed in the present paper. 
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